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SUMMARY 



The present paper describes the results of a stabil- 
ity investigation on symmetrical velocity profiles in a 
channel and of boundary- lay er profiles* The limitation to 
these two most important types of profiles was, however, 
not due to any 1 imi tat ion of our mathematical method . The 
effect of the friction was assumed to he vanishing and did 
not occur in the stability consideration so far as it had 
not been resorted to for preparatory asymptotic consider- 
ations. Proceeding on very general premises as regards 
the form of the velocity distribution, a proof was deduced 
of the elementary theorem that velocity profiles with in- 
flection points are unstable. Aside from this comprehen- 
sive theorem, there were obtained formulas of general va- 
lidity for the investigat ed . t ypes of profiles which dis- 
closed certain information relating to wave length, wave 
velocity, and amplification of the dangerous disturbances. 

Dynamically, the prof il es . with inflection points are 
identified, according to boundary- layer theory, by the 
existence of decelerating pressure gradients in such flows 
profiles, to which this particular stability investigation 
is inapplicable, may be encountered in slightly divergent 
channel s or on the hydro dynamical rear surface ( behind the 
pressure minimum) of cylindrical bodies, particularly if 
the body is sho rt . They may al so oc cur through sup erpo s i- 
tion of the velocity distribution without inflection point 
with a q.uasi- stat ionary vortex pattern, which can happen 
in the first stages of the formation of turbulence. The 
proved effect is , undoubtedly , of fundamental, importance 
for the creation of turbulence, 

* u Ein allgemeines Kriterium der Instabilitat laminarer Ge- 
schwindigkeit svert ei lung en . " Nachricht en von der Gesell- 
schaft der Wi ssenschaf t en zu G-ottingen (Mathematik) , .vol . 
I, no* 5, 1935, pp, 79-114, 
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IliTEODU'CTiOF '' * 

The t ransi tioxx of " lamina!^ flow"! with its clean, strat- 
ified layers of flow tubes to. strongly. intermingled, irreg- 
ular turbulent' "flow; constitutes one of th'e'most pressing 
problems of modern hydrodynamics- . It . is certain that this 
fundamental change in type of motion of the fluid is at- 
tributable to an instability in the laminar flow, for lam- 
inar flows of themselves would always constitute possible 
solutions of the hydrodynamib equat ions . The mathematical 
derivation of the expected instability of laminar flow is, 
as is known, beset with almo st insurmountable difficulties. 
While it is true th$,t the recent advances made on the sub- 
ject are very promising, it is equally true that almost 
every stability investigation made heretofore , relates to 
special laminar velocity distributions. There seems to be 
a palpable lack of 'general theorems in this field which 
afford ready classification of laminar velocity profiles 
according to their stability. 

The general instability criterion established herein- 
after, discloses a frequent and powerful mechanism of for- 
mation of turbulence. 

V 

The analysis is to proceed on the basis of two-dimen- 
sional velocity profiles, the velocities in this fundamen- 
tal flow- themselves to be very simply ' distributed; that is , 
t-o be largely dependent only on the coordinate transverse 
to the direction of flow. On this flow there are then su- 
perposed disturbances' which, according to the method of 
small oscillations,, are considered as small waves advanc- 
ing in the direction of flow. . The disturbances themselves 
shall also be two-dimensional*, while the' effect of fric- 
tion on the disturbances will, be disregarded to , a certain 
degree. The problem posed is, When can disturbances with 
increasing amplitud'e be superposed on the basic flow? 

This roughly outlined problem is by no means new. 
Back in 1880, Lord Rayleigh published an investigation 



*According to a note by H. B, Squire (Proceedings of the 
Royal Society of London A, 142, p. 621, 1933), the inves- 
tigation of the stability of flows depending only on one 
coordinate may, if the disturbance is three-dimensional, 
be simply reduced to the case of ' two-dimens ional disturb- 
ances . • • 
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which above everything else gave a re&dily obtainable nec* 
essarjr condition for instability. Subsequently, he again 
and again reverted to this same problem (reference l), and 
even other investigators, appreciating the importance of 
his findings* applied themselves to the problem, but with- 
out much success beyond the initial advance, No sufficient 
condition for instability has been set forth thus far. 



I. THE DISTURBANCE EQUATION 



let U represent the speed of the laminar fundamen- 
tal flow parallel to coordinate x and dependent only on 
the coordinate y at right angles to x. Being assumed 
two-dimensional, the disturbance can be derived from a 
stream function \|r (x.y.t), t being the time. The veloc- 
ity component of the disturbance in x direction is 3^/3y, 

3^ 

in y direction, - Thus the equation for the stream 

function of the disturbance inclusive of the friction 
terms - considering that U itself complies with the hy~ 
drodynamic equations and with limitation to linear terms 
in ijr, is: 

|^ +u m„ v A A t (1) 

at dx dy 3 3x s ' 

a ?\ z a 2 

where A is the Laplace operator + and V is the 

3 3y* 

kinematic viscosity. The equation is the vortex equation 
for the disturbance, which follows by elimination of the 
pressure from the hydrodynamic equations. A^ is the neg- 
ative vortex density of the disturbance, because A^ = 



3y " 3x ; 3^~ Sives the local time rate of change, U <|~r 



the. transport .of the disturbance turbulence by means of 

in 



the fundamental motion, - |^ the vortex transport 

the fundamental flow caused by the disturbance; V A A 
corresponds to the diffusion of the turbulence, by friction. 
The convective change in vortex density through the dis- 
turbance motion itself is neglected as nonlinear in. view 
of the as.sumed smallnesS of the disturbances. - 



Owing to the linearity of (l), the formation of q> 
from partial oscillations is permissible, which may be ex- 
pressed as 
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<p(y) e 



i(ax-pt) 



(2) 



that is, as waves traveling in direction x. In this equa- 
tion the complex method of writing is employed; for phys- 
ical application,' the real part of the above expression is 
to be selected. a is a real constant and equal to 2rr/X , 
where \ is the wave length of the partial disturbance; 
p may be complex equivalent to p r + ipi- The real part 
P r gives the cycle frequency of the disturbance, and the 

imaginary part (3^ the logarithmic increment , becoming 
positive for amplified, and negative, for damped- o sc illa- 
tions. The real part c r of the quantity c = | is the 
wave velocity. Writing (2) in (l) gives the disturbance 
equation for the considered partial o scillat ion: 



or, in nondimensional form (by measuring velocities in 
terms of a characteristic velocity of the basic profile, 
say, the maximum U ma2C , the length in terms of a charac- 
teristic width b of the profile), as 

(U-c) (cp»~cc s cp) - U" cp - - -~ (cp ,,u - 2a 2 cp»+a 4 cp) (4) 

CCR 

Eere we took the liberty of employing the same notation 
for the nondimensional quantities U, c, cc, etc. as before. 

U max b 

R is the Reynolds Number — . 

Slow suppose the kinematic viscosity is very small or, 
in more general terms, that the Reynolds Number R is very 
large. Then the experiment is suggested of suppressing the 
friction terms altogether to obtain an insight into the be- 
havior of the disturbance through a discussion of the fric- 
tionless disturbance equation 



and this is what Rayleigh actually did. He stepped from 
(4) to (5) without further discussion which, however, is 
naturally necessary. 




(3) 



(U - c) (cp" - cc s cp) - U» cp = 0 



(5) 
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II. ANALYTICAL PROPERTIES 03* . THE INTEGRAL OF THE 
DISTURBANCE EQUATION AND EEEECT OE INTERNAL FRICTION 



Starting with the analytical behavior of the solu- 
tions of (5), it is necessary to refer to aprevious re- 
port hy the author (reference 2). We assume c to "be 
purely real - that is, consider neutral oscillations, nei- 
ther amplified nor .damped. Now it happens, as we shall 
see later, that at one point within the fluid the there- 
existing fundamental velocity U "becomes equal to the 
wave velocity c of the considered partial oscillation. 
In other words, at this point a fluid particle always os- 
cillates in the same -d i sturhance phase. This point, called 
the critical point, is designated with subscript o. Now 
the origin of the y coordinate is placed at this criti- 
cal point, and the direction of' y is so defined that in 
the vicinity U - c > 0 for y > 0; that is, 

which is not to vanish, becomes positive. U 
pahle of expansion in power series around the 
point to: 



so that U 1 0 

is to he ca- 
critical 



U = c + 



y + 



2 



(6) 



The critical point U 0 = c represents a singular 
point of (5). Then two linear, independent solutions of 
(5) in the vicinity of the critical point can he repre- 
sented in the following manner: With P x (y) and P 2 (y) 

as power series in y* A and a as constants determined 
from (5), we have a fundamental system of solutions through 



epi = y Pi (y) = 



2U« 0 



y 2 + 



(7) 



CPs = Pg(y) + A qp x log y '= 1 + 



ay 2 + 



+ y log y 



(8) 



Hereby it is to be established once for all that 
log y for positive real y is purely real. When U" 
and likewise <p, does not disappear in the critical point, 
i.e., the transverse component of the disturbance, there 
results a- logarithmically infinite (disturbance) velocity 
in the x-direct ion, which is contrary to the smallness of 
the disturbance velo city assumed in the method of small 
oscillations. Regardless of the magnitude of the Reynolds 
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Numb erv the fr.iction mu'st ' "be taken 1 into ;aQcqunt . at .the 
critical point. The effect of the friction occurring in 
a small strip- (of; the order 'o;f magnitude of k (a R TJ 1 d ).~ 1/3 ) 
consists, first, in a flattening of the x component of the 
disturbance at the critical point - naturally it then re- 
mains finite - and secondly, in the appearance 'of a phase 
discont limit y in the disturbances . . Of the po ssible 
"branches of jfche log in (8), which would be mathematically 
available for "the continuation of the solution toward neg- 
ative y (naturally outside "of -the minute friction layer'),, 
the writer has, on the basis of previous calculations 
( reference 2) , chosen one which by means of the friction, 
is physically possible, namely, the analyt ical continua- 
tion of the logarithm by way of positive y through the 
lower complex - semiplane to negative y. When <p is rep- 
resented by ' 

U» 

1 ... + : T 7f- <P X log y (9a) 
o 

becomes 

l + ~ q> x (log |y| - irr) (9b) 

when y is negative. 

This transition substitution, however, signifies a 
phase discontinuity because, if the real part of 
i ( ccxi ' 13 1 ) 

9(y) e is selected for representing the disturb- 

ances, the. x component of the disturbance becomes 

/ U» \ 

- log y ...j cos (ax - pt) (10a) 

for positive y in the vicinity of the critical point, 
whereas ~ 

log |y| cos (ax - pt) + ^ tt sin (ax-pt) (10b) 

for negative y. The resiilt is a phase discontinuity which 
does not disappear even for the limiting case aE = 

In the present case we assumed c. purely real, .and 
only a very small amplification (Pi > 0) » that is, a 
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small imaginary part in c = & = c~ +- i ca . The zero 

a 

po iht o"f y "is placed a V the po int U =r~c , in accordance 
with the previous definition.. Attributing complex values 
to y and defining U from the power series (6), we again 
seek the singular point of (5) where the analytically con- 
tinued U ~ c = c r + i c i( This critical point U = c 
thus lies in the complex y plane slightly above the real 
, i c i \ 

axis (by — % — -I so that the real axis then has no longer 
a singular point of (5). Then at a sufficiently large 

/C| -l/3\ 

Reynolds Number Uji™ >(a R U'q) J, no friction effect 

on the real axis needs to be considered. Now we introduce 
an n coordinate shifted with respect to- y along the 
imaginary axis by placing the origin of r| in the critical 
point U = c r + i Ci which, for sufficiently small c^ 
gives 



y = r\ + ~i . (11) 



(See fig. 1.) 



Denoting the value of U f - and U" at U = c with 
U f c and U" c , the solutions cp x and Cp 3 in this case 

are written as before: 

<Pi = *n + 2t77 ^ + • • dza) 

u c . " 

cp s - 1 + . .. + & <p x ( T1 ) log r\ (I2.b) 

c 

We only use the values of the singular solution 
along the real axis y. For positive y, which is as- 
sumed to be large compared to the small quantity c^/u 1 ^ 
we have : 

U" 

q> = 1 + ... + =j2- <p!(y) log y (13a) 

o 

while for negative y, it' is: 

U» 



cp 2 = 1 + ... + - r S-cp x (y) (log |y| - ijr) (13b) 
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considering that, .the path* along which the logarithm must 
be analytically continued, lies "below the singular point; 
that is, the same transitional substitution -as for purely 
real c, which was' obtained through a boundary transit ion 
.from smair friction.* This is of particular interest to us 
because the neutral oscillation's are primarily considered 
as a limiting case of the amplified oscillations. 

Tor small damping, on the other hand, the transition 
is. from ; , 

U» 

?a = 1 + + ui7 cp i: (y) l0g 7 

for y positive to 

•• -. : .. ; "■ 

q). 3 = 1 + ... + --- qp x (y) (log.jyj + iii) 

' ' P 

for y negative.. 



(14a). 
(14b) 



Following the exclusion of . the vicinity of the crit- 
ical point through special- considerations, we have in <p 

and cp_ , the solutions of the f ri ct lossless disturbance 

s - * , • • 

equat ion , part icular solut ions of the complete disturbance 

equation for -very large a H.* 

The only remaining difficulty lies in the lowering of 
the order of the differential equation by 2 during the 
change from the complete to the f r.ict ionless disturbance 
equation. The result wasv of course., only two particular 
solutions and consequently noncompliance to all boundary 
conditions. If the fluid passes between two walls, for 
example, the tangent ial and normal component of the dis- 
turbance velocity must disappear at both walls. Limited 
to the frict ionless disturbance equat ion ( inclusive of the 
necessary correction at the critical point) , the stipula- 
tion is, of course, confined to the disappearance of the 
normal component only. But, according to the findings of 
various investigators (reference 3) the drop of the. tan- 
gential component of the disturbance to zero on the wall 
at high Reynolds Numbers actually takes place in an ox- 



♦Even in the case that Ci/U f 0 , while positive, is no 

— 1/3 

longer large compared to (a R U f 0 ) . the boundary 

transition to very largo R affords the same result when 
the writqr 1 s analysis (reference 2, p. 27) is repeated for 
this' case' " ' * ■ 



N.A.C.A. Technical Memorandum No. 792 



9 



tremely thin layer with the result that the effect of the 
wall ~f r-ict&on. become,s .consistently smaller • 

Besides this pair of asymptotic solut ions deduced 
from the frictionless disturbance equation, there exist 
rapidly decreasing ""boundary-layer-like solutions," the 
inclusion- of which in the frictionless solutions is neces- 
sary to insure complete compliance with the "boundary con- 
ditions. We shall give only the most elementary example 
of such a boundary-layer-like solution. Designating the 
distance from the wall toward the inside of the fluid y w 

and considering the particular case of c having positive 
real and imaginary part and he ing suf f ici ently great , 

, -1/3 



( 



_C 



^> la E 13\\ , IP W = value of U 1 at wall y w 



w 

0 j, the "boundary-layer-like solution, designated cp 3 then 
"becomes: 

^V^Rcy w 
9 3 - e 4 (15) 

The add it ion of this inwardly rapidly decreasing solu- 
tion necessary for compliance with the "boundary condition 
for the tangential component has but a minor effect on the 
distribution of the amplitude of the disturbance and the 
parameters a and c, as is readily proved. By suppress- 
ing the explicit occurrence of the internal friction in 
first asymptotic approximation for very large Reynolds 
Number s , the d et erminat ion of a critical R is foregone 
in favor of the desired generalization of the results. On 
the other hand, our investigation affords some particular 
information about the type of dangerous disturbances and 
constitutes for this reason a "useful preliminary f or .the 
disturbance problem with explicitly occurring internal 
f r ict ion« 

III. FORMULATION OP CHARACTERISTIC VALUE PROBLEM 
FOR SYMMETRICAL VELOCITY -PROFILES 



Following these preparations., we can f inally formu- 
late our problem* For tlie sake of simplicity, we first 
assume two-dimensional flow through channels as the fun- 
damental flow. The particular velocity profiles shall, as 
stated before , change very little in the direct ion of mo- 
tion; that is, U = U(y) , and the walls of the channel 
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shall he approximately parallel . Furthermore, the velocity- 
profiles shall he symmetrical with respect to the channel 
axis. At the wall, of course, U = 0 . Otherwise* the form 
of the" prof il es - is very little restricted. Profiles with 
an inflection point in-each half, stich as may occur in 
slight ly* divergent channels (fig* 2), are also included, 
hut ' not ' prof iles for which U itself changes sign. This 
includes profiles with return flow as well as separation 

profiles as a "limiting case (~ =0 at the wall\ At the 

• ■ Vdy / 

inflection point , of the included pro files (subscript s) , 

let U» s > 0, U"V< 0, so. that - U" is positive hetween 

the wall and point of inflection, and negative 'tetween the 
inflection point and the center. 

For these fundamental . velo citie s , the following bound- 
ary-value problem is to be solved: 

cp" - a 2 9 - cp = o (16) 

U — c 

with the boundary conditions that the normal . component of 
the disturbance, that' is, . Cp, " disappear s at both walls. 
The transitional' substitution suffered by cp at a singular 
point U = c in the case of neutral oscillations, has been 
previously established by a transition to the limit of van- 
ishing frict ion or disappearing amplification. The charac- 
teristic-value problem here is. rather unusual. The real 
parameter a is assumed predetermined; the quest is for c 
which occurs' in. nonlinear manner" in the differential equa- 
tion' (16). Problems of this' kind have boon so little ex- 
plored mathematically,, that the solution in question cannot 
be based upon general' existence theorems* The main task 
will be to establish necessary, and at the same time suffi- 
cient, conditions for U(y) , in order that complex charac- 
teristic values c may exist, for it is readily seen that, 
with a complex c as -character 1st ic value the conjugate 
complex c together with the conjugate complex character- 
istic function also represents a solution. Our interest 
centers on the solutions with amplification, for which c 
has a positive imaginary part. 

Since U is symmetrical, cp may be readily separated 
into a' symmetrical and an ant isymme t rical part which, indi- 
vidually, nus-t satisfy the . disturbance equation. As a re- 
sult , we need consider. onlj r half of the y-zohe, because 
then we have, either cp 1 = 0 or cp. =. 0 in the channel cen- 
ter in addition to the boundary condition cp = 0 at one 
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wall, depending on whether the distribution of. the disturb- 
ance amplitude . cp _ about , th.e central ;axis of the channel 
is symmetrical or ant i symmet rical •. This simplifies matters 
quite considerably because the problem narrows down to one 
singular point' in a semichannel* 

The .particularly important .:. pro blem o.f* stability of 
boundary— layer, profiles, is. treated* in a subsequent section, 

: IV.. EAYLEI&H 1 S EQUATIONS 



His principal result,., summed up briefly., is a necessary 
condition for the occurrence of- amplified oscillations. 

Writing the differential term on the left-hand side of 

(16) 

<p» - a 2 cp - ----- cp = L(cp) (17) 

U — c 

we form the integral term 

2b 

/ [<P L(q>) - 9 L(9)] <ly w (18) 

o 

between the boundaries of, the y^zone, the two channel walls. 
2b =. channel breadth. • Conjugate complex quantities are 
overlined. Therefore, 

1(9) = <P" - a 2 cp - SE1L. 9 (17a) 

U - c 

whence Cp and cp shall vanish, conforming to the boundary 
co nd it ions at the zone boundaries • The problem is to find 
when a nonvanishing imaginary part of c is possible. On 
this premise there is no singular point on the real y axis 
so that the integral summarily gives: 

2 b 

- 31 oi / Tz-S^-pr 9 ? ^w < 19 ) 

o p.~ c i 

This term then must vanish for solutions of the differen- 
tial equation besides the boundary conditions. If U u 
does not change • sign ,. this is » however , ' not po ssible for a 
0^ other than zero , so that amplified o scillat ions are im- 
po s sible except for pro files with ■ an inflect ion' point • The 
condition, however, being only' necessary/ Rayleigh 1 s theo- 
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rem merely states that profiles without an inflection point 
reveal no instability in the sense used here.* Any state- 
ment as to the "behavior of profiles with an inflection 
point is as yet impo ssible . 

Another fact "brought out by Lord Eayleigh is that, by 
the eventually jneutral oscillations (c^ = 0 , c = c r ) , the 
wave velocity c p of the partial oscillation must be equal 
to the basic velocity at one point, thus ever assuring .the 
existence of a point with TJ - c = 0 within the fluid. We 
adduce a proof for this fact , which is somewhat more simple 
than Rayleigh's. (See reference 1, vol. VI, 1913, p. 199; 
also vol. I.) 

Let us assume that c > U max (maximum value of U) , 

so that no singular point exists within the fluid. We com- 
pare 

q>" = (a 2 + g-?~)<P ( 1 6) 

with 

f - -JEL. f (20) 

that is, compare two solutions of (16) and (20) which dis- 
appear at the wall (y w = 0) . In addition, let Cp 1 = fi 

for y^ = 0. Multiplication of (16) by f and (20) by cp 
followed by subtraction of the results and integration 
starting at one wall gives: 

f <p ! - ft cp = / a 2 f <p dy w (Si) 
o 

So long as <p and f are positive, it follows: 

cp ! 4? i 

..9 > ~ and (p > f (22) 

after which, the solution of (20) becomes: 

^w d v« 

f = (c - U) / (23) 

o (U - c) 



*Many profiles without inflection points may, by applying 
a correction explicitly containing the internal friction - 
that is, in the second asymptotic approximation for ver?r 
large Reynolds Numbers, become unstable, as proved by cal- 
culation on special velocity distributions. The amplifi- 
cation as a result of the small friction correction is, of 
course, of lower order of magnitude than that found here. 
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f is therefore always positive and greater than zero on 
the other wall,; consequently, cp cannot disappear on the 
other wall, according to (22) or, in other words , it is 
necessary that c < U ma x # * 

By the same argument, c must he greater than the min- 
imum value of the fundamental velocity (zero hereafter). 
The proof is the same, word for word, if we take 

y"w dV«r 

( »- •>/ r^V 

for f. It should he noted that for the proof of Rayleigh 1 s 
two formulas, the assumption of symmetry of U was not 
ut il iz ed ♦ 

Rayleigh, "being unahle to advance beyond these two 
necessary conditions for amplified or neutral oscillations, 
decided to approximate the steadily curved velocity profile 
by a profile consisting of straight pieces (polygonal pro- 
file), in order to side-step the mathematical difficulties. 
The transitional conditions to be fulfilled by cp at the 
bends of the profile were easily established. They corre- 
spond to the condition of equal no rmal component of the ve- 
locities and equal pressure on both sides of the bend. So 
when Rayleigh** divided the profile into three strips of 
constant vortex density (d^/dy = constant), for example, 
he arrived at an ordinary quadratic equation for c and 
was able to show that in case of re-entrant corners, ampli- 
fied oscillations are always present. But the rough approx- 
imation of the velocity profile through a polygon suppresses 
essential parts of the process. The profiles treated here- 
inafter are always curved. 

V. EXISTENCE OE NEUTRAL CHARACTERISTIC OSCILLATIONS 

We begin with the neutral oscillations which are possi- 
ble in the case of a symmetrical channel profile, and which, 
according to Rayleigh, insure a singular point c = 0 with- 

*That c < Uj^x for the special case of U ,! <0 (pro- 
files without an inflection point), is readily seen by 
forming the integral term 

2b _ 

/ [cp L(cp) + cp L(cp>] dy w . 

o 

**See reference 1 (Vol. Ill, p. 17). 0. Tietjen's discus- 
sion on boundary-layer profiles is similar. 
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in the/ zone. If • cp is other than zero in this critical 
point « v say,- equal to Cp 0 - the previously cited transi- 
tional substitution (9) and (13), upon change from positive 
to negative y, that is, on passage from channel center to- 
ward the wall,., suddenly creates at Cp 1 an additive compo- 
U» • 

nent ' - — ~ i rr cp 0 , whereas <p remains continuous. It will 

U o ■■ 

be shown that such an abrupt increase is impossible with a 
neutral characteristic Oscillation. To this end, we consid- 
er the expression: 

'' ■ n & - *r £ 

Incidentally we add that (24) is proportional to the 
time average of the product of both disturbance velocities; 
that is, it is intimately connected with the momentum trans- 
port of the disturbance. Since cp in the case of neutral 
oscillations satisfies a differential equation with real 
coefficients, both cp r and cp ^ must individually sat isfy 
the differential equation. Consequently, (24) is constant 
so long as no singular points occur in the dif f er ential 
equation. As a result of the cited sudden growth of 9 
on passing through the singular point U = c, (24) also in- 
creases suddenly upon transition from greater to smaller y 
through the singular point, and by an amount 

U" 

as is readily seen from a simple calculation of the behavior 
of Cp and cp 1 at the singular point. Since, according to 
the boundary conditions, the term (24) vanishes at both the 
center and the wall, thero is no possibility for such a sud- 
den increase. Consequently, either cp must be equal to 0 
at the singular point for the neutral characteristic oscil- 
lations, i.e., only the regular solution <P X supplies the 
neutral characteristic oscillation, or else the singular 
point must drop out altogether, whence U n 0 = 0; that is, 
the phase velocity must equal the fundamental velocity TJ g 
at the inflection point, provided, of course, that the 
profile has one. 

Starting with the possibility of cp = 0 at the criti- 
cal point, we compare , * 
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with 



9 n = (a s + —J--) q> (16) 

*v-.(r?t) f ' • (20) 



and compare those solutions <P X of (l6) with the solutions 
f x of (20) , both of which disappear at the' critical point 
and whose power expansion 'otalrts at the critical point with 
y. Then , 

Multiplying (16) by f x , and (20) by cp x followed by 
integrating from the critical point .( subscript o) , gives 

y 

f, i^i - fi qp'i = " / a s f x cp x dy - (26) 

o ■ 

f x ..and cpx are, in accord with the assumption for nega- 
tive y (from critical point toward 'the wall) negative at 
first, and the term on the ' right-hand side , pp s it ive • So 
long as f x and cp x are negative, therefore, 

f f , CD 1 

_i Jk_ and y 1 < f ( 27 ) 

x i SP X 

Since f x is still negative at the wall ^= - u?~~) 

rather than' passing through zero, cp x can so' much less 
disappear at the wall unless ■ ' c' = 0. ' Then, of. course, hoth 
f x and <p x disappear at the wall .at which; the critical 
point- then falls. 

Ap.plyi.ng the above argument In. the case c = 0 to the 
zone between wall (y w = 0) and center (y w = b) , gives . 
for this zone, • 

v • 

9i f i ~ qpi f ! i = / Ot' 2 - f x 9i iyw > 0 . ■ (28) 



.0 



' " q>i > f x - 5 f- ( 29 ) 

. . ■ : . u :w 

so long, as . f x = -2- =" and . cp x . are positive. Since 

u o ^ w 



16 ff.iA.O.A. Technieaa Memo randum. JJp ♦..79-2 

in, the center fi = nr~-" 0 ' i;fc follows from ( 28 ) 

,.JL_ cp[ x >o, that is, q> ! x > 0 in the center. According- 
ly w 

ly, no cp ? x corresponding" to" a 2 > 0 can satisfy the 

"boundary condition at the center. This condition thus com- 
pels a = 0 at the center which, in conjunction with the 
previous condition c = 0 leads to the regular solution 

f ± = ---- or with other notation, to cp = 'U, which' remains 
IT w 

as the only one for a characteristic solution and actually 
is such a one, according to the properties of U. 

As the only characteristic oscillation existing in the 
same manner for profiles with or without an inflection 
point, we obtained the abnormal oscillation -with a = 0, 
c = 0;' that is, with infinite wave" length and vanishing 
wave velocity and with the symmetrical amplitude distribu- 
tion cp = U. This solution itself has been known for a 
long time} we prove here its singularity. 

However,, fundamental velocity profiles with an inflec- 
tion point may also have a . neutral characteristic oscilla- 
tion for which c =U S (fundamental velocity at the in- 
flection point). Mathematically, the question is as to the 
proof that the Sturm- Lionville formula 

qp" + x 9 - T 7— V~ <p = o 

u - u s 

with the .conditions Cp = 0 for . y w = 0 and Cp f = 0 for 
y w = b,- • aside from the. infinitely many positive character- 
istic X values exist ing according to general theorems, 
has also a negative characteristic value \ = - a ? (fig. 3). 

The proof proceeds -from the solution cp = U ± U s for 
a = 0, which has disappearing tangent in the center.- This 
boundary condition in the center is maintained for the so- 
lutions studied at a 4= 0 . The addition of a 2 > 0 slows 
up the drop of these solutions from the center toward the 
wall, with the result that the nodal point Cp = 0 of these 
solutions, which for a = 0 lies at the inflection point 
of the fundamental velocity, finally approaches the wall; 
that is, we have the desired characteristic oscillation 
(fig. 4). The conclusion of the steady .dependence of nodal 
point position oh the parameter (& s ) Is, as known,. made 
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in proof of the oscillation theorem (reference 4). Since 
the parameter (here a 2 ) may hay.©,.- only . positive values , it 
summarily follows that always one, and only" one , such sym- 
metrical disturbance flow is possible, "because characteris- 
tic oscillations with several nodes cp = 0 can exist only 
for a negative parameter. The existence of the characteris- 
tic solution could, moreover, also be proved by repeated 
application of our previous deduction. We call this charac- 
teristic solution cp s and put its parameters- a = a s , c = U s « 

Any doubt as to the possible occurrence of an antisym- 
metrical disturbance flow for c = U s with 9 = 0 in the 
center, is refuted in a subsequent section. 

Lord Rayleigh himself had pointed to another possible 
solution for U" = 0,' According to the frictionless dis- 
turbance equation (5), it might be conjectured that at such 
a point a solution U = c existed; that is, a disturbance 
confined only to the critieal layer U = c itself.' Consid- 
ering the high values of the differential quotients of <p 
in such an oscillation, the friction must in any event be 
taken into consideration. Since L. Hopf (reference 3) has 
proved that these oscillations are damped, they are merely 
mentioned, 

VI. EXISTENCE OF AMPLIFIED CHARACTERISTIC OSCILLATIONS 



Their existence under certain conditions will be proved 
by considering solutions adjoining the previously deduced 
neutral characteristic solutions. Our method is vaguely 
reminiscent of the well-known perturbation theory of the 
characteristic values (reference 5) with, however, substan- 
tial modifications imposed by the prevalent singularities. 
Rather than proceeding, as customary, from the known to the 
hypothetical solution, we attempt to double back from the 
hypothetical (amplified) solution to the known (neutral) 
char a c teristic solution. 

Let cp n represent the neutral characteristic solution 
and c n , a n its parameter s , the parameters c ( complex) 
and a being assumed. One solution (pj is to satisfy the 
disturbance equation with these parameters without, however, 
being a characteristic solution, whereas q)j is to comply 
with a boundary condition, say, that on the wall. In this 
manner cpj is defined up to one factor; c, a, and cpj are 
subsequently appropriately determined as being adjacent to 
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c n , a A5 and cp n . ■ .. . 

Subtracting. ..the dif f erent ial ..equat ion for . 'Cp n ; . •• . 

(U, - c n ) (9 n n ~ CC 2 n cp n ) - U" cp n ,= .0 ' (-30 ) 

from that for <Pj: 

(U- c) (qp«j - a 2 qpj) U" cpj = 0 (31) 

•give's . 

(U - c) - cp n )» - a 3 (cpj - <p n )} - -U» -(.qpj - <p n ) 

■ l - . y ■ (32) 

= (c - c n .)cp», n + |(U - c) a 2 - (U - c n ).a 2 n j q> n . . 

(<Pl - 9 n )" - a 2 (cpx - cp n ) - _#U""<<Pr - 9n) = g (33) 

with 

g== oTtf - 9 n " U~- c °^ n 9x1 * a 9x1 (34) 

g contains . only .cp n ,. and the disturbance parameters* Hav- 
ing assumed c as complex, we have in reality no zero 
point of U - c , that is, no singular point in the differ- 
ential equation for .cpj - cp n . . Be it- noted that th/e coeffi- 
cients of the.- homogeneous part cpj - <p n - of. (33:)- correspond 
to those of the differential equation for <Pj. 

Then with <$n as a linear solution of 

9" - a s cp ^ T T"— <P = 0 (31) 

independent of Cpj, we can write cpj - <p n in form of 

Vi - <P a = qpj / W S <Pu dy w - cp ri / W .g qpj dy w + C cpj"' (35) 
o • o 

with integration constant C in which, hy normalizing 

: Vx V'li - cp^ cp I]: = - i 

The integral is to extend from the wall as subscript w 
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indicates. The problem then reduces to finding the Neces- 
sary and sufficient condition which satisfies the "boundary 
condition at the cent er. through r cpj , that is , cp' "j =='■ 0 or 
else Cp'j - cp f n =■ 0. With b as the distance of the wall 
from the channel center, we have for y w = "b, according. t;o 
(35): 

b b 
cpVj, . cp' n ^ tp-Vj / g,^ f , dy w cp'-n / dy w + Cqp'l 

0 0 

.:■ ;• : ■ ■ - V . ■ • (36) 

To deduce a necessary condition, it is assumed that 
■''■I." 0 in the center',. -so» that .".'q> r j j becomes-, other than 
zero in view of the presumed linear independence of cpj . 
Then, of course, it follows' that 

b 

/ e q>j. dy w = 0 (37) 
o 

in accordance to (36). This go nd it i on , suf f ici ent since 
for y w = b, according to' (36) 

" b 

' <P'l. I 1 ~ C te / . g ; 9ll ^w] = 0 (38) 
o 

as well as 

b 

q>i [1 - C ~ J g <Pii dy w ]' = cp n : (39) 

" ■ ■ o : 

according to (35). 

But cp n is not zero in the center because there 
cp* n = 0, whence the term in the brackets must likewise be 
other than zero. Consequently, . . <p* j.= 0 in (38). 

.• How we prove the existence- of solutions with ampl.ifica- 
tionfor profiles with inflection points by means of (37), 
for it permits the calculat ion : of c" (at: least near c n ) 
as function of a (near a n ) . ■ 

•We first., consider the vicinity . of cp n = U. f c n = a n - 0* 
The solution cpj adjacent to this neutral characteristic, 
solution- of. the disturbance equation with small parameters 
c (as sumed complex- and, with pp sit iye' imaginary part ) and. 
a (real) is established through . qpj = 0 at the wall and' 

the further condition that the derivative' of cpj at the 
wall equals that of Cp^ . that is, , cpVj. == XJ 1 ^ for. y w = .0 . 
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Iia this- case 



g = ■■ e „.u» +■ a 3 U • " (40) 
u — c 



and 

g qpj = c cp» x 4- a 2 (U - c) q> z (41) 

by having recourse to the differential equation for cpj. 
Accordingly, the condition 

/[c cp n ! + a 2 (U - c) cpj] dy w = 0 (42) 

o 

must be met for a characteristic solution, according to 
(37) .* 

The evaluation of (42) proceeds on the derivations of 
several approximate representations cpj • ' It does not suf~. 
fice to approach q)j through cp n , because the operation 

must be effected in the dangerous proximity of a singulari- 
ty. We again expand according to r\ from 'the singular 
point L T = c = c r + i cj, so that r\ is approximately 

equal to y - ^?JL. As "before, subscript c denotes the 

values of U ! , etc. at the singular point; c r is chosen 

positive. The expansion is then made for a range extending 
on the one hand, from the singular point to the wall, where 

c U", 



n approximates — a — **— , and on the other hand, to 

1 ~~ U« c 2U* C 3 

a small positive y = e ; the latter (c) to be chosen so as 

°r 

to make c 2 small relative to — r» whereas G is large 

c t .. ! L cj 

relative to yyj — ;— f Here as well as later, it is simply as~ 

I u c! 

sumed that c^ proves of smaller order of magnitude** than 
c r . Accordingly, we assume G as being about equal to 

c 3/4: 

-F— 5 C ? X X is to equal U ! w at the wall. The expan- 



*This condition (42) can also be proved direct with Cp f j = 
0 in the center by substituting ( U- c) cp^j-U" Cpj for 
a 2 (U-c)cpl according to (31). 

* *Thi s relat ion ship between c r and c\ is already indicated 
by the fact that the roughest approximation for Cpj, namely, 
II , gives an expression for c r a c cording to the integral 
equation (42), but not for c^. 
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sion, starting from the singular point, gives 



U" rt c U" c 2 _ U" f c " o 
u 'w = V*a ™r ° 2 + V c 

The development of cpj with two unknown coefficients A and 
B is: 



By determining A and B from the initial conditions, we ob- 
tain for fPj in the cited zone up t o terms of the approxi- 
mate order of c r 2 : 



r U" / C \ >j / U" \ 

»i ={i - « 57V [!•«(- — ) + i]}(Ti'cn+ -5»tf) 

c c 

+ C c ~ ° s u^O C 1 + ^ loe ^ (43) 

The imaginary part of log ( -* ] is i tt , while 

v u c / c* 

that of log Tj at y = C is equal to i — which, ac- 

cording to the definition of € is very small. 

The region of y = € up to the channel center y w = h 
requires a representation of the maximum imaginary part in 

c • 

cpj . Since c ^ 7fM*"T"' we can for this range, develop 

u c 



according to Cj_, which affords as a first approximation: 



-ar x c^ + cc 



U - c r (U - c r ) 



= 0 



iSTow it is possible. to set up a fundamental system of. 
solutions, one of which satisfies the initial condition 
cp = 1 , 4> f = 0 for y = e. the other cp = 0, cp 1 = 1. In 
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the absence of singularities the development may "be effect- 
ed according to parameters i c^ and a 2 , and it is found 

that in this fundamental system the first imaginary compo- 
nents go with ci* The indicated development of cp starts 
with 

<P = Qx(jr) + i Ci Q 2 (y) + a 2 Q 3 (y) 

where the Q has real, restricted values which are not de- 
pendent on parameters Ci and a 2 . Q 1 may he explicitly 
expressed, according to the early statements about (20) . 
By suitably combining these fundamental solutions, the 
joining to (43) may be effected for y = €. Here again, 
Cj_ is to be of lower order of magnitude than c r » Remem- 
bering the remark about (43), according to which the imagi- 
nary part of log f- ~— J equals - rr i , while that of log 
r\ disappears for^ y^= c e, the maximum imaginary part* in 
9j is easily obtained as 

U" 

c r Tr-S- i it U( y ) (44) 



TP 



c 



How' (42) can be evaluated: qpj changes for- c = a = 0 
steadily to U. With this rough approximation (42) gives 
the equation for c r to the first approximation: 

b 

U? w + c r ^'i^) + cx 2 / U 2 dy w = 0 



c. 



If, however, the condition of orthogonality is com- 
plied with, then <p f i(b) = 0, whence 

c r = -f- J U 3 dy w (45) 
u w o 

that is, actually positive, which' proves the assumption 
about c r as being in accord with the condition of ortho- 
gonality. To define Cj_ demands closer approximations 



♦The statement about tyj can also be confirmed insofar as 
an explicit expression for cpj can be given, at least for 

a = 0, that is, qpj. = - U f w c(U-c) J as is read- 

o (U - c) 

ily seen when written in the differential equation* 
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(43) or (44) for . q>j. There is* obtained 

' " yn" b ' r - : U" 

: oi - a 2 c r 77 ^3 it" / tJ 2 - dy w = c r 2 77 V tt (46) 

Equations (45) and (46) contain' the necessary and at the 
same time sufficient conditions for the : existence of ampli- 
fied char act eristic solut ions-. It remains to be established 
whether these conditions conflict with the assumptions rela- 
tive to the disturbance parameter made incidental to cpj # 
As to (45) , this has been proved; as to (46) , only the fact 
that c^ > 0 was utilized in the representation of cpj 
should be noted. Equation (46) can be complied with only 
when the profiles have an inflection point, where U !, w > 0, 
For the others, U ,! w < 0, so that we obtain a contradic- 
tion in accordance with the first formula of Rayleigh. And, 
as the assumptions for c r and c^ are also confirmed, 
the existence of ampl if ied char act eristic oscillations in 
the case of profiles with inflection points is proved. It 
is seen that c^ increases go much more with c r as, un- 
der otherwise identical conditions, U n is greater at the 
wall. 

Although this very fact proves the instability of such 
profiles, it should be very instructive to construct a 
neighboring solution with amplification for the neutral 
characteristic solution cp s with the parameters c n = U s 
and cc n = cc s . Again denoting this particular solution' by 
Cpj and the parameters near U s and a Q by c, a, we put 

c - U s = Ac = Ac r + i c?i 

(47) 

a s - a s s = Aa 2 

Ci again being assumed positive, so that (37) suitably 

transformed with the aid of the differential equation for 
9 S > becomes : 

b b 
/ ——4 <L T f—~ cp s CP! dy w + f La? cp s cpi dy w = 0 (48) 
o (U - c) (U - U s ) o 

Now q> s is visualized as normalized, so that it is ex- 
actly equal to 1 at the inflection point. If y is the co- 
ordinate starting at this point, the development of cp s 
starts with 1 + ky where k is constant, which need not 
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"be known further. We. may regard <Pj a.s . det ermined "by the 
fact that it "becomes 1 at the singular point U c = c = 



U s + A c r + 



i c. 



and satisfies the "boundary condition 



<Pj = 0. We again effect the development about the singular 



point in terms of r\ - 7 

interval - € ^ y ^ + C, 
obtainable "by choosing £ ~ 

■ The development of cpj 



_ Ac 

where 
Ac 



that is, at first, for an 



AC 



U' 



3/4 

u r c' 

according to 



Tmt e 3 <€ 



Ac_ 



U», 



T| starts with 



(49) 



while 



c. 



Cp s = 1 + k 



Ac 



+ k T) 



Owing to the s t eady ' t ran s it idn from' CPj into 
ishing Ac and Aa 2 , the constants k and k x , 
by the condition at the wall, 'differ "but little. 



(50) 



for van- 
det ermined 
Then the 

integration, according to (48) for the interval of y be- 
tween ~ € and + €, gives for the first integral in (48) 
(limited to terms at least linear in Ac): 



. c 



i TT 



(51) 



or, when considering that 
the inflection point y w = 



U H c through' development from 
s becomes e qua 1 to ^ M s ^ c , 



U' 



s 



where U in c 



i s assumed negat ive and U r s po sit ive . 
U l!t 



_ „.__s__ c ..i JL + 
u> s s 



U"'s V 
77% c 



Out side of this int erval <Pj 



r rr i (52) 
may again be developed , 



C is, on the one hand, chosen so small that a few terms 



of a development of 9j 
|yj < g, that |U - U s 



according to r\ 
^>Ac for ly| 



suffices for 
> C and u 



U - c 



and Cpj can be developed according to the parameter Ac. 
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according to the parameters Ac and Aa s • The first ap- 
proximat ion ( linear in Ac), .for the missing constituent of 
the first integral of (48) is obtained by substituting Cp- 
for cpj in this zone: 

L o (U - U s ) s+€ (U - U s ) J 

From the determination of e , it follows that the tran- 
sition e — >0 can he effected (up to terms which, like e, 
disappear) within the "brackets without changing the value. 
We put : 

E=lim{ / € n —r 9s 2 dy w + / _El _ (p s a dy w } (54) 
C==0 1 o (U-TI S ) s+G (U-U s ) ^ 

E contains only known quant i t ies and remains f init e, because 
the constituents which become infinite, cancel out. Conse- 
quent ly , 



U»i b , yitt ^ 

E Ac r - — -I" CiTT+Aa 2 / cp s 3 dy w +i|ciE+ — TrAc r |=0 (55) 
^ s ^ XJ g 



The explicit interpretation of this relation requires 

the knowledge of E and / <p s 2 dy w , which would call for 

the calculation of the proven neutral characteristic solu- 
tion cp g . But the existence of a neighboring characteris- 
tic oscillation is readily apparent without it. To illus- 
trate, assume E < 0 as can be proved (section VIIl)« Then 
the imaginary part of (55) gives: 

U ,,f o A c r , e „ x 

ci = - ~ s -2 - tt X 56 ) 

U' S S E 

Since c ^ according to assumption is to be positiv^ 
and U tn s negative, A c r should in this case, be negatives*. 
Writing c^ according to (56) in the real part of (55) 
gives: 



which proves that a satisfactory solution of the necessary 
and sufficient condition for amplified characteristic o s- 
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cillations is possible as soon as Aa is chosen negative* 
It can also he proved that for E - 0, in which A c r - 0 
or E > 0, an amplified characteristic oscillation would 
exist as soon as A a 2 is negative. Continuation of the 
characteristic solutions over a in the zone around a g 
is therefore only possible in one direction, namely, for 
decreasing a, where E should also have the same sign. 
The order of magnitude of c^, according to the above for- 
mulas is at least as great as that of A c r ^' The result is 
the start of the development of function c(oc) at two 
points, namely, for a = 0 and a = a s in general formu- 
las. 

VII. EQUILIBRIUM OF BOUNDARY- LAYER PROFILES 

A particularly important case is now analyzed as to 
stability. Idealized to a certain extent to suit our pur- 
poses, the boundary-layer profiles are to be so defined that 
the velocity U rises from 0 at the wall to U maX at 

distance 6, and then remains constant to infinity (fig. 5), 
At y w ~ 5 > a discontinuity in TJ" may be permitted', for, 
integration of the differential equation for Cp over a 
small interval around y w = 8 , shows that nevertheless, in 
addition to cp , cp* also must remain continuous at this 
po int . 

Now a fundamental sy st em of cp solut ions can be set 
up for the zone U = constant = U maX , namely, cp = e" r -^ 
and cp = e +a 2w, in which, when a is positive, the second 
may not occur because of its infinite growth in. a charac- 
teristic oscillation. The characteristic function in the 
zone from y w = 8 to y w = oo behaves like e ~ a3r w. Our 

analysis shall be restricted to the zone 0 = y w = 6 , so as 

to riile out the infinitely remote point, which is an essen- 
tial singular point of the differential equation. Then the 
charac t er ist ic solut ion, aside from cp = 0 at the wall 
(yw ~ 0) , must satisfy the other boundary condition 

cp r + a cp = 0 for y w = 6 (58) 

in order to make the joining with the cited solution form 
e - f X7 w s p 0 ssible . 

The neutral charac t eristic solut ions are e stabli shed 
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very similarly to those for profiles in the channel, and 
jfrotn an analogous conclusion it follows that 0 =5. c = U max 
and that as fat as profiles without an Inflect 1 off" point 
are concerned, ohly the neutral characteristic solution 
<p=U, c = a = 0 exisis* ftd construction proceeds from 
the solution of the differential equation for c - U s , a = 
0 , which sat isf ies the boundary condit ion f or y w = 8 • 
But this solution itself is simply U - TJ S . Changing to • 
a 3 > 0 while preserving the value of cp at y w = 5 "by 

normalizing, the q> value is increased twofold: first , as 
an increase of cp for y w ~ 6 corresponding to the "bound- 
ary condition cp 1 = - a cp upon advance to smaller y w ; sec- 
ond, as effect of a on the coefficients of the differen- 
tial equation, as known from the channel profiles. Finally, 
when a becomes large enough, cp ' ceases to he negative at 
the wall (fig* 6), For this a = a s , where q> = 0 at the 
wall, the characteristic solution cp s exists. 

How we construct their neighbor ing amplified character** 
istic solutions. The previous formulas for the channel pro- 
file are in any case in the neighborhood of <p n = U, c = 
a = 0 not to be transferred by simple limit process to the 
boundary-layer profiles, because there it gave, for in- 
stanc e , 

2 h 

c r = T ~ ~ / U 2 dy w (45) 
u w o 

so that for boundary-layer profiles, where b goes to in- 
finity, the integral would become infinite. For this rea- 
son, a new integral condition is derived for the existence 
of an adjacent characteristic solution, the analysis again 
being restricted to the finite zone 0 < y % 6 * 

Subtracting the differential equation for the neutral 
characteristic solution cp n from that for the adjacent so- 
lution Cp j , we obtain 

(<Pj - cp n )» - cfCep! - <p n ) - (q>i - <p n ) = 6 (^3) 

with 

« = %~~r ^ " ihHr a3 » ^ + a * ^ ' (34) 

It is s impl est to derive a neces sary condit ion that <p be 
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a characteristic solution., as, follows: : The differential 
equat ion f o r cp j • . 

is written symbolically L(cf)j) = 0, so that the above equa- 
tion for cp j — cp n becomes: 

L(cpj - cp n ) = g 

■ 

Next we formulate: 

6 ' - ' >, 

/ |<Pt L(q>! - cp n .) - (<p z - <p n ) L(qp z ) | dy w 




(59) 



•which with consideration to (58) gives the necessary condi- 
tion: 

5 

/ <pi g *yw = '(<% - Cq>i q>„)_ _ s (so) 

That this condition is also sufficient, is seen in similar 
fashion as in section VI. We form 

cp f ! - qp' n + <x(q>l - 9 n ) 

S 8 
= (9' I fa cp x ) / g.cp n dy w -(9 f il + a ^il) / * ?! ay w +C(<p' i+a <Pj) 

0 o 

(61) 

for Cp j - 9 n at y w = 8 according to (35). 

Considering the. condition (60) looked upon as fulfilled, 
-and Cpjj expressed as 

cpj cp f IX - cp T ! cpjj = - 1 

the right-hand side becomes: 

5 

(cp 1 ! + a cpj) / g cp n dy w + 0(cp f I + a <Pt) - (a - a n ) <p n 
o 

+ (<P ! ! + a cp x ) ( a - a n ) <p IX <p n 
and the whole equation (61) may be written as: 
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f 6 ^ 

(9 f i+acpi) -jl-C - / g <¥>n dy w +(cc n -cc) cp lT q> h j - 9 f n +a n cp n =0 (62) 

The multiplier of qp'j ■ + a <pj on the left side is differ- 
ent from 0, "because according to (35) and.. (60), it is: 

( 5 ' ' > . 

<pj |l - 0 - / g cp I3[ dy w + (a n - a) q> n <P n j » .q> n (53) 

for = 8, and <P n (8) * 0; otherwise <P ? n(S) == 0 for 

y w = 8 » that is, <p n ^ would have to disappear,* " Conse- 
quently, the "bracketed term does not disappear and 
9*1 + a Cpj = 0 for y w = 8 , according to (62). 

The determination of cpj near ' cp n .= U, c = a = .0 as 
hefore., through = o and cp 1 j = U T W - at the wall, is 

■followed by the ' e valuat ion of the newintegral condition' 
(60) : 

8 ' 

/ ^c cpnj + a s (U - c) q> x Uy w = (a n - a) (cpj <P n ) . 

O ^ «> W~^ 



as 



- c r U' w = - a U max 2 



or rr TT 2 

a u max 



c r = ~;f*±- ( 54) 



w 



A comparison of this formula with the previous one re- 
veals another power of a at the 1 imi t ing • t rans it ion from 
channel to boundary-layer profiles. The amplification is 
(deductions "as- "before, especially with formula (44)): 



a U" c 2 U» 

Ci —g. - TT U max = — ~g- (o5) 

u w u w 



The relationship existing between both formulas for 
c r (45) and (64) becomes so much clearer when comparing 
symmetrical channel profiles, which are to have a variable 
U from y w = 0 to y w = 8 , but constant tJ = Umax over 
a great er 1 ength 2b x around the c enter , with the boundary- 
layer profiles (compare third profile in fig. 2). Cpj is 
for the channel profile in the* vicinity of constant TJ ap~ 
proximately: 

♦That cp n (8) •* 0 is also readily apparent from the cited- 
construction of the neutral- characteristic solutions. 
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e a(bi+6) e -ay w + e -a(bi+8) e +ay w 
u max " ^7^™-^^ 
b 

The constituent / a 2 U 9j dy w - in (42) becomes equal 
to a U max 2 tanh a b x ^ and approaches a 2 U max 2 bi for 
small a b x according to (45) , and approaches a U max 2 
for large a hi in accordance with (64). The validity of 
(45) is therefore contingent upon a b x <€ 1, which for 
boundary-layer profiles, where *b x : — ^oo is impossible, no 
matter how small a may be , ■ in which case formula (64) • 
holds for c r . 

As to the adjacent amplified oscillation for 9 S> 0L = 
a s , c = U s , suffice it to state that (56) and (57) are 
preserved; only the there-existing integral needs to be ex- 
tended to b = 

VIII. ILLUSTRATIVE EXAMPLES 



These examples treat the neutral characteristic solu«* 
tions cp s , previously referred to, as well as the formulas 
for wave velocity and amplification in the vicinity of a = 
0 and a = a s . First we assume sinusoidal channel pro- 
files. Thus, 

. U = U s + (U max - U s ) sin (65) 

The velocity U s at the inflection point is then connected 
with s, the distance of the inflection point from the 
wall , through the relation: 



which follows from U = 0 for y =0, whence the differ- 

w 

ential equation for a characteristic solution at c = U g 
rtins as follows: 

1 TT 2 

q>» ~ a 2 <p + 7- rsr 7 cp = 0 (68) 

{ b - s) 4 

To allow at the same time for the boundary condition 
at the wall, we put Cp ==.sin E?W; p is determined from 

. b ' 
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tlie "boundary condition for y w = h. Fo.r a, the differen- 
tial equation gives: 



a 3 = 4r 



( 



1 17 



(69) 



Since 0 H s < g and-' a 2 is always positive, p may as- 
sume only the value tt/2 in the center, according to the 
"boundary condition cp* = 0, tut not the tt value, corre- 
sponding to 9 = 0. The result is therefore- the character- 
istic function 



• TT 

8 Si 



(70) 



and 



b 2 • 4 



C 1 - 0" 



5- - U 



(71) 



for s = ^, ah = 0.559 tt , for example (fig. 7). 

o 

For the amplified adjacent solutions with this pro- 
file in the vicinity of 

c r = 0,462 a 2 h 2 U. 



= 0 , a = 0 ; 



max 



Ci = 7.59 



U. 



c r 


M 


r max 




c r 


= u s = 


c r 


= 0.061 


c i 


= 0.194 



5/2 



(72) 



U 



i- = 0.336 - 0.82.1 



(73) 



max 



^max 

where ah must he smaller than 0.559 tt. 

We proceed now to the sinusoidal houndary-lay er pro- 
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"files, for which 

U - U s + (U max - U s ) sin ^ for y w i 5 

6 " S 2 (74) 

U = U ma£ for y w > 6 

The differential equation for cp in the zone 0 ^ 

y ^ 6 runs as follows: 
w • 

cp" - a s 9 + 3- — <P = 0 (75) 

(S - s) 4 

with the "boundary conditions: 

cp = 0 for y w = 0 

cp* + acp = 0 for y w = 8 

a was assumed positive. By making 9 = sin — the 
other "boundary condition gives: 

a = - p cot p (76) 
The differential equation furnishes the further relation: 

_L_ 1 = -JE (77) 

s 2 sin p 

5 

We consider only positive values of p, as .negative p 
would give no new characteristic function. In view of 
equation (77), sin p > 0. Therefore, 

10 ' z p 



sin p 



and, since 



tl g — ±. — n. < n 

2 i _ s 2 * 



p must he smaller than - n, according to (77). As a is 
to he positive, p = tt/2 according to (76). The solution 
of the transcendental equation (77) for several spec if ic 
cases , follows: 
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For s = 0., p = 3j-, .. a = 0 
for' s = |, p = 0.658 tt .(= -^|~ ij). - 0.356 n 

"for s = | , p = 0.737 tt (= ~f |q ~ tt), .. a6 = 0.675 tt 
The characteristic function "becomes:- . 

. py OT 



(78) 



<J>< 



sin — : 



. PS 

sin — — 
6 



for 0 y w 8 



sin p a§ -ay w , > . 

e e for jr - 6 

ps w - 

sin r 



(79) 



For a closer descript ion of the shape of ' Cp we 
. » ... s 

point to the position of the maximum at ?-w = JL. Figure 

' 5 2p 

8 shows - plotted for s == A. 



For' the boundary-layer profile with s = ^ near 



= 5. 

3 



c r - 0 , a = 0, . the amplified neighboring solutions are: 



c r - 1.025 a 5 U max 

max . max x max ' 



(80) 



in the vicinity of ' c r =' U a = 0.414 'U ma2c , Ofi = 0.356 tt 
. . ( c r : = 0.083 a* 5 2 Unax + 0.312 U max 

c ± = 0.250 U max - 0.604 o" r (81) 

6^ c„ . 

■ = 0.2 80 - 0. 67 5 



'max u max 
where a8 must he smaller than 0.356 tt. 
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Figures 9 and 10 show c r versus a and ver sus 

c r according to (72), (73), (80), and (81), The start 
and finish of the curves are shown as heavy lines connect- 
ed by dashes after interpolation, soas to give an idea of 
the probable amplification. 2Io conclusions are drawn 
therefrom. If it is preferred to obtain the entire distri- 
bution of the cited curves the attempt at generally appli- 
cable formulas must be foregone in favor of a calculation 
fitting the part icular • case • . And that itself is easy, 
once the difficult problem of the existence has been ex- 
plained. 

Translation by J. Vanier, 
National Advisory Committee 
for Aeronautics 
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Figure 2. 



Figure 4. 
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